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A SIMPLE METHOD FOR ESTIMATING TERMINAL VELOCITY 


INCLUDING EFFECT OF COMPRESSIBILITY ON DRAG 


By Ralph P. Bielat 


SUMMARY 


A generalized drag curve that provides an estimate 
for the drag rise due to compressibility has been obtained 
from an analysis of wind-tunnel data of several airfoils, 
fuselages, nacelles, and windshields at speeds up to and 
above the wing critical speed. The airfoils analyzed had 
little or no sweepbaek and effective aspect ratios above 6.3* 
A chart based on the generalized drag curve is presented, 
from which the terminal velocity of a conventional airplane 
that employs a wing of moderate aspect ratio and very 
little sweepbaek in a vertical dive may be rapidly esti- 
mated. In order to use the chart, the only data that 
need be known about the airplane are a low- speed drag 
coefficient, the wing critical speed, and the wing loading. 
The terminal velocities for three airplanes were computed 
in order to Illustrate the use of the method and chart. 

Good agreement between the estimated terminal velocity and 
the measured flight terminal velocity was indicated for 
all three airplanes . 


INTRODUCTION 


Several high-speed military airplanes in dives have 
encountered difficulties that could not be easily con- 
trolled by normal means. These difficulties, which may 
consist of diving moments, large changes in trim, large 
stick forces, tail buffeting, and the like, occur in high- 
speed dives when the speed of the airplane exceeds the 
critical speed by a large amount. For those airplanes for 
which maximum diving speeds are at or near the critical 
speed, little or no- trouble occurs. The more recent 
fighter airplanes, however, have terminal Mach numbers 
well in excess of the critical Mach number and, as a 
result, often encounter difficulties in dives. Determi- 
nation of the terminal velocity of the airplane is there- 
fore Important' in order tha t the. T pr obabi li ty of encoun- 
tering trouble In d i ve s ,)0 m£?y t ima t e d .. 
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The terminal velocity is ai30 important because it 
forms the outer limits of the V-G diagram. Usually the 
outer limit of the V-G diagram is established by multi- 
plying the maximum level-flight speed of an, airplane by 
an arbitrary factor somewhat greater than 1.0. The termi- 
nal velocity of most recent airplanes, however,- generally 
falls much below this arbitrary maximum speed, and these 
airplanes are therefore unnecessarily penalised by extra 
weight because they are designed for conditions that are 
not reached in actual flight. 

The present report outlines a simple method for 
obtaining the terminal velocity of an airplane in a verti- 
cal dive and includes an estimate for uhe drag increase 
due to compressibility effects. The wind-tunnel test data 
were obtained from model tests conducted in the Langley 
2k- inch and 8-foot high-speed tunnels. All the data 
presented herein were obtained for zero lift. 

The problem of determining the terminal velocity for 
airplanes for which the terminal velocity is near the 
critical speed is comparatively simple inasmuch as a 
constant value of drag coefficient can be assumed. The 
diving speeds of most present-day airplanes, however, 
occur beyond the critical speed and the problem is nert so 
simple. The following two factors are involved* (1) the 
determination of the critical speed and (2) the rate- -of ' 
drag increase at speeds above the critical speed. 

The critical speed used herein was arbitrarily ’taken 
as the critical speed of the wing-root section. Pressure- 
distribution data obtained from wind-tunnel tests -were 
used to determine the critical speed, which is defined, as 
the flight speed at which sonic velocity is reached 
locally. If experimental data are not available however ,- 
the methods outlined in references 1 and 2 can be used 
for the determination of the critical speed. Selection, _ 
of the critical speed at the wing-root section for use in 
terminal- velocity estimation is justified on the. gremmds. 
that the root section usually has a lower critical speed 
than any other component part of the airplane. 'The wing 
root has the lowest critical speed because of . .its., high 
thickness ratio and contributes a large part of the total 
airplane drag because of the large fraction of the wing, 
area concentrated at the inboard sections of tap-ered, wings . 

The rate of drag increase at speeds above the., criti- 
cal speed is more difficult to determine in the calculation 
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of the terminal velocity than the critical speed. A 
study of the drag of airfoils, fuselages, nacelles, and 
windshields ha3 been made from wind-tunnel test data in 
order to determine the effects of compressibility on the 
drag. Eecause the rate of drag increase at speeds above 
the critical speed is so great, it was found that, within 
the accuracy required for terminal- velocity calculations, 
an average rate of drag increase may be used. A curve 
indicating the average rate of drag increase is presented 
herein. This curve was derived from an analysis of wind- 
tunnel data. 

The method described herein for obtaining the termi- 
nal velocity of an airplane in a vertical dive has been 
in use at the MCA since 19I4.I. Publication of the method, 
however, had been delayed pending the investigation of 
constriction corrections to be applied to ihe wind-tunnel 
data and the completion of high-speed dive, tests made 
with several airplanes in order to compare terminal 
velocities obtained in flight with terminal velocities 
estimated by the simple method described herein. This 
method is not applicable to airplanes that utilize wing 
shapes of low aspect ratio and large sweepback but should 
be applied only to airplanes of conventional design that 
employ wing shapes of moderate aspect ratio and small 
amounts of sweepback due to wing taper ratio. 

SYMBOLS 

V velocity 

a speed of sound in air 

M Mach number (V/a) 

Cj) drag coefficient 

C L lift coefficient 

p mass density of air 

S wing area 

W weight of airplane 

p atmospheric pressure at any altitude 


k 


FAC A ACR No . 


Y ratio of specific heats (l.ipO for air) 

t/c ratio of thickness to chord of wing 

Sub a cr i pt 3 s 

cr critical (when local sonic velocity has been 

reached on some point of body) 

min minimum 

T terminal 


DESCRIPTION OP NOBELS 


Airfoil models .- The airfoil models used herein 
represent two classes of airfoils - namely, the conven- 
tional NACA sections and the more recent low-drag high- 
critical-speed NaGA sections. The conventional NACA air- 
foil sections are characterized by pressure distributions 
that have high peak pressures occurring near the leading 
edge. The low-drag NACA airfoil sections have compara- 
tively flat pressure distributions with the peak pressures 
occurring at approximately 60 percent of the chord behind 
the leading edge . 


Airfoils typical of the conventional airfoils are 
the NACA OOO9-65, 0012 - 65 , 23012 , and 23OII4..7 sections; 
a current transport-model airfoil that has an NACA 2215 
section at the root and tapers to an NACA 2212 section at 
the tip; and the Davis airfoil with a thickness ratio of 
20.15 percent. The low-drag airfoils include the following 
NACA airfoil sections: 


16-215 

16-509 

16-515 

47-215 
65(218) -220 


^ 5 -type modified, — = O.I96 

66,1-115 
67- lib.. 5 
67,0-215 


The effective aspect ratio of the airfoil models 
tested varied from 6,5 to infinity. 


Fuselage models.- The fuselage models are typical of 
fuselage shapes in use on current airplanes. The various 
fuselages represent bomber, fighter, and transport air- 
planes. Figure 1 shows the side-view drawing and the 
fineness ratio in side elevation of the different fuselage 
shapes. These fuselage models were tested in conjunction 
with wings (shown as das h ed li nes in fig. 1) and represent 
a wide variation in ftfgf interference. 
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Nacelle and windshield models .- The data for the 
various nacelles and windshields we re obtained from refer 
ences 3 and I;., respectively. The nacelle and windshield 
designations used herein correspond to the designations 
used in references 3 and I 4 .. All the nacelle models were 
tested with the same wing model, which consisted of the 
outboard panel of a wing section designed for use on a 
bomber airplane. The wing was a thick .low-drag airfoil 
that had an MCA 65(218) -221 section at the root and 
tapered to an WACA 06 ( 2x15 ) -lj.l 6 section at the tip. The 
windshields were tested with a wing-fuselage combination. 
Drawings of the nacelle and windshield models are shown 
in figures 2 and 3 * respectively. 


RESULTS AND DISCUSSION 
Drag Characteristics 


Drag analysis .- In order to obtain a correlation of 
the rate of drag increase at speeds above the critical 
speed, the drag results for the various component parts 
of the airplane have been reduced to nondimensional param- 
eters: that is, ^ D / C U m in P^o^ed a S a i^st M/M c ~. 

for each part tested. The use of these parameters 
represents a convenient method of making the data non- 
dimensional in such a manner that the unknown quantities 
are expressed in terms of the known quantities. 

The drag results at speeds up to and above the criti- 
cal speed for the conventional MCA airfoils are presented 
in figures l\. and 5* Figures 6 to 9 show the variation of 
°D/°D mm Klth M / M cr for the low-drag high- critical- 

speed airfoils. It will be noted that all the airfoils 
presented in figures 4 , 6 , 7 > and. 8 exhibited approxi- 
mately the same rate of drag increase at speeds above the 
critical speed; for this reason a curve of the average 
rate of drag increase at speeds above the critical speed maybe 
used. An average increase in drag of approximately 30 per- 

cent above the minimum drag was indicated at ~ — — 1 . 0 ; 

M cr - 

at speeds of only 10 to 15 percent above the critical 
speed, however, the drag increased. approximately 90 
to 200 percent. This rapid increase in drag at speeds 
above the critical speed is associated with the formation 
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of compression shock waves and their effect on the ■boundary- 
layer over the surface of the airfoils. The family of 
airfoils used in figure 5 showed less percentage of 
increase in drag at the critical speed than the NACA 0009, 
0012, or the low-drag high-critical-speed airfoil sections. 
Both published (reference 5) and unpublished high-speed 
data show that the NACA 230- series airfoils differ from 
most of the other airfoils in that the critical speed can 
be exceeded by as much as 0.15 in Mach number before any 
serious changes in the aerodynamic characteristics of the 
airfoil occur. The critical speed of the NaCA 230-series 
airfoils is therefore exceeded by approximately 7~ percent 

£L. 

before the same percentage of increase in drag occurs as 
is shown for the other airfoils. The importance of this 
difference in the shape of the drag curve above the criti- 
cal speed on the estimation of terminal velocity is dis- 
cussed in the section entitled ''Terminal-Velocity Calcula- 
tion." 

The rapid increase in drag before the critical speed 
is reached, which is shown for the NACA 67-.Hl>.. 5 airfoil 
in figure % is due to early separation of the flow over 
the after portion of the airfoil. This condition also 
affects the method for estimating the terminal velocity. 

An error in estimating the terminal velocity when the 
flow separates will occur only for those airplanes for 
which terminal velocities are at or near tne critical 
speed; this separation of flow will not appreciably affect 
the determination of the terminal velocity for high- 
performance airplanes for which the terminal velocity 
occurs at speeds well above the critical speed. 

Figure 10 shows the variation of with 

M/M cr for several fuselage shapes and fineness ratios. 

The drag increments for the nacelles and windshields are 
presented in figures 11 and 12, respectively. The criti- 
cal speeds for these bodies were based on the wings with 
which the models were tested and were determined for the 
wing-root juncture. The effect of compressibility on the 
rate of drag increase at speeds above the wing critical 
speed for these bodies is similar to that for the air- 
foils , 

In the correlation of the average drag increases of 
the various components of the airplane throughout the 
Mach number range, a generalized drag curve was derived 
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and is presented in figure 13 . The data presented in 
figures [j. , 6, 7? Si 10, 11, and 12 were used to obtain 
the generalized drag curve. The generalized drag curve 
is an average of the drag data for the airfoils, fuselages, 
nacelles, and windshields at speeds up to 10 percent 
above the critical speed. Only the average drag of the 
airfoils at speeds from 10 to 15 percent above the criti- 
cal speed was used. The generalised drag curve was extra- 
polated by use of a straight-line extrapolation from 15 
to 25 percent above the critical speed. The straight- 
line extrapolation is believed, to be sufficiently accurate 
for estimation of the terminal velocity in this region 
'where the drag rises rapidly due to compressibility 
effects. 

C onstriction corrections .- Corrections for constriction 
effects have been applied to the data. The constriction 
corrections have been determined from pressure measure- 
ments obtained in the Langley Zip-inch and 8-foot high- 
speed tunnels on NACA 0012 airfoil models of various 
sizes. The magnitude of the corrections, applied to the 
drag coefficients amounted to less than one-half of 1 per- 
cent of the dynamic pressure q at low speeds and 
increased to approximately 2 percent of q at the criti- 
cal speeds and to approximately 5 percent of q at a 
value of the Mach number below the choking speed of the 
tunnel. The corrections to the Mach numbers amounted to 
approximately one-half of these values. The constriction 
corrections were such that the coefficients were reduced 
and the Mach numbers were increased by the values stated. 
The greatest percentage of increase in correction, as 
would be expected, occurred for the models that had the 
largest ratio of model area to tunnel area. 

Comparison with flight data .- Figure lip shows the 
variation of over-all drag coefficient with Mach number 
for the XP -51 airplane as measured in flight and the 
variation with Mach number of the wing-profile drag at 
the mid-semi span station measured by the wake-survey 
method. These flight data are preliminary as corrections 
to the data have not been applied. The results obtained 
by use of the generalized drag curve in estimating the 
drag increases with Mach number are also shown in fig- 
ure lip for comparison with the flight measurements of 
over-all drag and wing-profile-drag data of the XP-51 air- 
plane. The curves for the wing-profile drag and the over- 
all airplane drag in flight begin to rise rather steeply 
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at about the same Mach number. This fact tends to justify 
the assumption that the wing-root critical speed is a 
suitable criterion to use in terminal-velocity calcula- 
tions. The estimated drag derived from the generalized 
drag relation indicates higher drag coefficients at Mach 
numbers of approximately 0.55 to 0.75 than are shown for 
both the measured wing-prof ile drag and the over-all drag 
coefficients. Of more importance, however, is the good 
agreement that is shown for the values obtained by use of 
the generalized drag curve and the measured flight data 
at Mach numbers greater than 0.75* which is the region 
where the terminal Mach number usually occurs. 

Figure 15 shows a comparison of measured flight drag 
and estimated drag for the XF2A-2 airplane of reference 6. 
An important difference in the drag curves occurs at Mach 
numbers around the critical Mach number. The estimated 
drag indicates lower drag coefficients than do the flight 
measurements. This difference is believed to be due to a 
combination of early shock formation on the cowling and 
airplane-wing roughness, which is believed to have caused 
some separation of the flow. Good agreement is indicated 
between the flight measurements and the estimated drag in 
the region where the drag coefficients rise steeply, which 
is the region that determines the terminal Mach number. 


Terminal-Velocity Calculation 

The generalized drag curve (fig. 15) may be used as 
an approximation in determining the terminal velocity of 
an airplane in a vertical dive. The terminal velocity is 
reached when the drag of the airplane is equal to the 
weight of the airplane . The drag of the airplane in a 
dive combines both airplane and propeller characteristics. 
In the present analysis, however, zero propeller thrust 
is assumed and the propeller drag or thrust is therefore 
neglected. At supercritical speeds the drag or thrust 
caused by the propeller is considered to be negligible as 
compared with the drag of the airplane, particularly if 
the pilot throttles the engine and adjusts the propeller 
to a high blade-angle position. 

The terminal velocity for zero-lift conditions is 
given by the relation 
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v _ FTT 

T p p'Od 


( 1 ) 


or, in terms of the terminal Mach number My with the 
speed, of sound equal to \l~^, 


Mm — 


V, 


T 


IW 2 1 _1_ 
Srp% 


( 2 ) 


Equation (2) can be rewritten in the form 


„ , W 2 1 1 

Mm ~ ~ ~~ ~~ 


S Y P %in Cd / C ^1 


n 


w 


1 


Y pSC 3min C D/°D m 


m 


(5) 


where 


psc n 


is constant for each airplane at the 


'min 


altitude for which Mm is calculated and 


_2a_ = /Vn 

"Alin 


which is obtainable from the generalized drag curve. 
Equation (3) can then be solved for the parameter 
W 

- 7 — . Figure 16 shows solutions of equation (3) for 

. sS0 %ln 

various values of M cr and 1 


pS %in 


In computations of the terminal velocity, the only 
data that must be known about the airplane are the minimum 
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drag coefficient at zero (or approximately zero) lift 
coefficient or a low-speed drag coefficient whereby the 
minimum drag coefficient can be computed by use of the 
generalized drag curve, the wing critical speed, and the 
wing loading.- Values of these quantities, all of which 
are used in calculations other than those for the terminal 
velocity, are easily obtained. With these values known 

W 

for a oarticular airnlane , the parameter can oe 

pSC %in 

calculated for different altitudes; then, for given values 


of M 


or 


and 


W 


pSCr^ 


the terminal Mach number can be 


mm 


obtained by use of figure l6. 


In order to illustrate the method of obtaining the 
terminal velocity graphically, the terminal velocities 
have been calculated for the XP2A-2, P-39E-1, and P— 2+7 air- 
planes. The pertinent data for these airplanes, are given 
in the following table: 


TABLE' I 
AIRPLANE DATA 


Airplane 

___ __J 

M or 

C D . 
min 

W 

S 

(lb/sq ft) 

XF2A-2 

0.6l (flight) 

.66 (corrected) 

0.022 (flight) 

26.1 

P-39N-1 

0.675 (estimated) 

0,018 (estimated) 

34.1 ! 

- - - j 

P-4? 

0 . 61+ (wind tunnel) 1 0.020 (flight) 
.6$ (corrected) j 

i 

45.0 ; 

1 

1 


By use of these data the parameter 


v* 


PS- c d, 


is 


mm 


computed for each airplane . The use of figure l6 to 
estimate the terminal Mach number is illustrated for the 
P-i+7 airplane at 15,000 feet altitude. The variation of 
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terminal Mach number with altitude thus obtained for the 
three airplanes is presented in. figure 17. 

Also included in figure 17, for comparison with the 
estimated variation of terminal Mach number with altitude, 
are records of flight data for the XF2A-2, the P-J4.7, 
the P-I4.7C-I-RE, and the P-39N-1 airplanes. The flight 
record for the P-I4.7C-I-RE airplane was obtained by the 
late Major Perry Ritchie in a terminal-veloci ty dive made 
at Wright Field in July 1943* The points represented by 
circles were obtained from a dive of a P-47 airplane made 
by a test pilot for the Republic Aviation Corporation. 
Unfortunately, a complete dive history is not available 
for this dive but it is believed that, had one been 
available, it would have followed a path similar to that 
obtained by the la,te Major Ritchie for the P-li?C-l~RE air- 
plane. It is further believed that the test points 
obtained at altitudes of 22,000 feet and 10,000 feet 
represent entry into and pull r out from the dive, respec- 
tively. Data for the XF2A-2 airplane were obtained from 
reference 6 and the data for the P-39N-1 were obtained 
from dive tests made at Ames Aeronautical Laboratory. The 
present method for estimating the terminal Mach number 
yields results that compare favorably with the flight 
measurements; the difference between the two is no greater 
than 0.02 in Mach number. This method for estimating the 
terminal Mach number is therefore believed to be suffi- 
ciently accurate for usual engineering purposes. 

The section entitled "Drag Characteristics" indicates 
that the NACA 230-series airfoils and airfoils similar to 
the NACA 230-series could exceed the critical speed by 
approximately 0.05 to 0.15 in Mach number before any 
important changes in the aerodynamic characteristics 

occurred. At ~4— = 1.0, therefore, the NACA 230-series 

^cr 

airfoils and similar airfoils did not show the same 
percentage increase in drag as was shown for almost all 
the other airfoils and for the generalized drag curve. 
Since in the calculation of the terminal velocity the 
critical speed of the airplane is based on the critical 
speed of the wing, it can be expected that for airplanes 
utilizing NACA 230-series airfoil sections or similar 
sections the estimation of the terminal velocity will be 
in error. If the generalized drag curve is used in the 
estimation of the terminal velocity, the Indicated wing 

-• •. 'JMK 1 

critical speed must be increased approximately percent 
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for the NACA 230 -series sections. This correction was 
applied to the critical speeds of the P-k7 and XF2A-2 air- 
planes (see table I), since these airplanes have NASA 230 - 
series sections. The dashed curve on figure 17 for 
M C r = 0 . 6 ij. is the result obtained if the indicated criti- 
cal Mach number is used rather than the effective critical 
Mach number, which is about 7— percent higher. 


Langley Memorial Aeronautical Laboratory 

national Advisory Committee for Aeronautics 
Langley Field, Va. 
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Elliptical cross section Circular cross Section 



Figure 2.- Nacelles tested. 
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Figure 4 .— Variation of ratio C D / Co m ; n with Ari/Ar} i 
tor NACA 0009 and 00/2 airfoils. 
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Figure. & Variation of ratio C D / C Dm - n with M/fv1 cr 
•for three Conventional MAC A air/bi/s . 



Figure 6 . — Variation of ratio ^o/C 0/riin with tslj/ui cr 
for NACA 16-509 and 16-5/5 airfoils. Transition 
fixed at 0./O chord . 





(p) Transition fixed at 0.10 chord- 
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(b) Transition fixed at 0.40 chord. . 

Figure. 7 Variation of ratio C D /c D . with 
for several NACA airfoils . 
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Figure & ■- Vanatior of ratio Cp/Cp^n with M/Mcr 
For several NACA airfoils . No trarsit/or) . 
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Fiqure 9.- Vernation of ratio ^ol^orpin with /\slj f\ 
for Davis airfoil and sev&rat NACA airfoih . 
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Figure // .“ Variation of ratio ^o/^D m , n with hsi/i\r cr 
for several nacelle shape s. 
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Figure 14 Variation o-F measured -flight drag arid 
estimated drag with hdach number For XP-51 
airplane . 



Mach number , A 4 

Figure 15.- Variation of measured flight drag and 
estimated drag with Mach number for XF2A-2. 
airplane . 
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Figure 17.- Comparison o-P estimated va riaiion and Plight' variation oP terrain a/ /vfach 
number with attitude for several airplanes . 
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